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ABSTRACT 

T-joint  configurations,  commonly  found  in  aerospace  applications,  provide  a 
challenge  for  structural  health  monitoring  (SHM)  strategies  that  employ  wave 
propagation  based  interrogation  techniques.  The  challenge  arises  from  the 
complexity  of  wave  propagation  through  the  T-joint  interface,  and  the  consequent 
mode  conversion  and  reflections.  This  work  presents  a  numerical  characterization 
of  undamaged  T-joint  using  simplified  two-dimensional  (2D)  plane  strain  finite 
element  (FE)  models,  where  excitation  from  a  piezoelectric  actuator  is  simulated  on 
one  side  of  the  T-joint.  The  incident  and  consequent  reflected  waveforms  are 
separated  in  the  frequency-wavenumber  domain,  and  spatially-integrated  multiple 
component  mode  coefficients  are  formulated  in  the  frequency-spatial  domain  with 
the  intent  to  characterize  wave  reflections  and  mode  conversions.  Mode  coefficients 
from  the  undamaged  configuration  are  compared  to  T-joints  with  blind  side  notch 
damage  located  in  the  fillet  radius.  Trends  from  the  2D  FE  numerical  studies  are 
verified  by  a  three  dimensional  (3D)  Scanning  Laser  Vibrometer  test  setup  for  the 
undamaged  case,  which  consists  of  an  aluminum  T-joint  for  selected  fillet  radius 
(0.6  cm)  and  stiffener  thickness  (0.3  cm).  The  setup  enables  3D  full  wavefield 
surface  measurements  of  the  front  and  back  sides  of  the  plate,  plate- stiffener  radius, 
and  stiffener.  Recommendations  of  how  to  use  the  mode  coefficient  formulation  as 
a  damage  quantification  tool  for  the  T-joint  are  given. 

INTRODUCTION 

The  stiffener-plate  configuration,  or  so-called  T-joint,  is  ubiquitous  in  aerospace 
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components,  and  most  commonly  represented  in  the  wing  airfoil-spar  and 
bulkhead-stiffener  interfaces.  A  nominal  fdlet  is  generally  designed  at  the  interface 
to  alleviate  potential  stress  concentrations,  which  can  produce  initiatory  cracks 
under  high  cycle  loading  conditions.  However,  these  T-joint  configurations  provide 
a  challenge  for  structural  health  monitoring  (SHM)  strategies  that  employ  wave 
propagation-based  interrogation  techniques,  such  that  partial  transmission, 
reflection,  and  mode  conversion  are  modified  by  the  degree  of  the  fillet  radius, 
excitation  frequency,  material  constituency,  and  stiffener  cross-sectional 
thickness[4].  Since  the  success  of  this  proposed  SHM  strategy  to  assess  damage 
severity  largely  depends  on  the  quantification  of  the  mode  conversion,  a 
fundamental  understanding  of  the  wave  interaction  due  to  the  T-joint  presence  is 
necessary. 

Previous  work  has  been  conducted  on  understanding  T-joint  configurations 
of  various  material  constituencies  and  a  range  of  damage  quantification  techniques. 
Modal  techniques  employed  by  Ooijevaar  et  al.  [1]  and  Ghoshal  et  al.  [2]  used 
scanning  laser  vibrometry  to  detect  damage  produced  in  composite  plate  with  two 
T-shaped  stiffener  sections  and  induced  long  saw  cuts  on  an  aluminum  wing  panel, 
respectively.  In  terms  of  wave  propagation  techniques,  Greve  et  al.  [3]  analyzed 
Lamb  wave  interactions  of  a  cracked  welded  joint,  and  rolled  and  plate  girders 
using  2D  and  3D  FEM.  Simulations  and  experiments  showed  cracks  that  are 
oriented  perpendicular  to  the  direction  of  wave  propagation,  cause  strong 
reflections.  Oppenheim  et  al.  [4]  modeled  a  2D  simulation  of  Lamb  wave 
transmission/reflection  at  a  T  joint  with  no  fillet  radius  using  Abaqus  and  3D  model 
of  plate  girder.  Experimental  measurements  using  PZT  sensors  located  on  the  web 
in  pulse-echo  mode,  and  strong  So  mode  reflections.  Similarly,  Morvan  et  al.  [5] 
measured  Lamb  waves  by  a  laser  vibrometer  in  a  T-joint  generated  by  a  contact 
piezo-composite  transducer,  and  normal  component  of  the  surface  wave 
displacement  of  the  plate  and  stiffener.  Sohn  et  al.  [6]  analyzed  mode  conversion  of 
T-joint  with  collocated  PZT  before  and  after  the  stiffener  using  a  2D  simulation 
under  plane  strain  assumptions.  Other  studies  concluded,  with  experimental 
evidence,  that  damage- sensitive  features  are  successfully  extracted  as  damage- 
induced  Lamb  wave  modes  have  unique  characteristics  compared  to  those  modes 
produced  by  the  stiffener  or  varying  thickness  [7,  8], 

The  objective  of  this  paper  is  to  further  understand  the  mode  conversion 
mechanism  at  the  T-joint  interface  using  frequency-wavenumber  filtering 
techniques  per  Ruzzene  [9],  both  as  a  function  of  the  fillet  radius,  with  and  without 
blind  side  damage,  and  with  experimental  validation  that  provides  high  fidelity 
spatial  resolution  using  3D  laser  vibrometry. 

MODE  COEFFICIENT  FORMULATION 

This  analysis  is  restricted  to  the  frequency  range  of  the  fundamental  symmetric  and 
antisymmetric  Lamb  waves,  So  and  Ao,  respectively.  As  illustrated  in  Figure  1,  it  is 
assumed  that  a  pure  So  or  Ao  mode  is  generated  to  interrogate  the  structure.  The 
interaction  of  an  incident  So  mode  with  a  stiffener  located  causes  a  partial  So 
reflection  and  partial  mode  conversion  into  the  Ao  mode,  while  simultaneously 
producing  a  partial  So  and  Ao  mode  transmission  into  the  stiffener  and  the  plate.  As 
previously  described  by  Ayers  et  al.  [10],  upon  filtering  the  data  in  the  frequency- 


wavenumber  domain  and  windowing  the  data  in  the  time-space  domain,  w  (x,  co)  is 
the  filtered  response  in  the  spatial-frequency  domain  of  a  single  extracted  mode  at 
frequency  to. 


SSA 


Figure  1.  Schematic  of  T-joint  and  mode  conversion  identification 


Laser  vibrometry  techniques  allow  for  various  points  to  be  scanned  over  the 
distance,  L.  From  the  measured  responses,  the  value  of  the  mode  coefficients 
may  be  estimated  over  a  spatial  domain.  This  allows  the  estimation  of  its  spatial- 
independence,  and  more  importantly  of  a  spatially  averaged  value,  which  is  less 
affected  by  amplitude  fluctuations  due  to  noise,  or  by  general  trends  associated 
with  geometrical  spreading.  In  compact  matrix  notation,  the  mode  coefficients 
may  be  expressed  by  the  following  for  an  interrogating  So  mode,  but  is  likewise 
easily  applied  for  an  interrogating  Ao  mode: 


n,(w)  = 


(1) 


where, 


(2) 


%jM  = 
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and, 
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where  R,  T,  and  C  parameters  represent  the  reflection,  transmission,  and 
converted  coefficients,  u  and  w  are  the  in-plane  and  out-of-plane  displacement 
components,  respectively,  and  x,  are  the  integration  spatial  coordinates.  It  is 
important  to  note  that  this  spatial  integration  method  departs  from  existing 
literature,  where  the  mode  amplitudes  are  taken  at  a  fixed  location  [11], 

APPLICATION  TO  SIMULATED  DATA 

The  mode  coefficient  formulation  is  applied  to  numerical  parameterization  for  the 
undamaged  and  damaged  configurations  using  simplified  2D  plane  strain,  where 
mode  conversion  estimation  is  shown  as  a  function  of  the  fillet  radius,  excitation 
frequency,  and  stiffener  cross-sectional  thickness  (Table  I  and  Figure  2). 


Table  I.  Test  Matrix  of  2D  FEM  Parameterization 


Set 

Signal 

Plate-Stiffener  Geometry 

Notch  Damage  Geometry 

Frequency 

Fillet  Radius 

Stiffener 

Thickness 

Depth 

Orientation 

[kHz] 

[mm] 

[mm] 

[mm] 

[deg] 

1 

200.0 

[0.0:3.0:12.0] 

3.0 

- 

- 

2 

200.0 

[0.0:3.0:12.0] 

3.0 

3.0 

0.0 

3 

200.0 

6.0 

3.0 

[0.0:2. 0:8.0] 

0.0 

4 

200.0 

6.0 

3.0 

3.0 

[0.0:15.0:75.0] 

I 

L. 


a)  (b) 

Figure  2:  Schematic  of  2D  FEM  T-joint:  (a)  Profile  with  identified  boundary 
conditions,  (b)  Magnified  view  of  FEM  parameters  at  fillet  with  notch  damage 

From  Figures  3(a-b),  the  mode  coefficients  from  a  no-damage  configuration  vary 
parabolically  as  a  function  of  the  fillet  radius.  With  a  constant  damage  depth  in 
Figure  4(a),  the  mode  coefficients  varies  linearly  as  a  function  of  the  fillet 
radius.  As  the  fillet  radius  increases  (greater  than  6  mm),  the  mode  coefficients 
from  undamaged  to  damaged  are  approximately  equivalent  (Figures  4(b)).  The 
thickness  of  the  fillet  radius  (greater  than  6  mm)  requires  that  the  damage  depth 
does  not  penetrate  the  plate  thickness,  and  sensitivity  is  negligible.  The  %sc-  mode 


coefficient  is  sensitive  to  the  increasing  damage  depth,  whereas  the  Qo  mode 
coefficient  shows  perturbation  (greater  than  4  mm)  once  the  damage  depth 
penetrates  the  plate  thickness  (Figures  5(a-b)).  The  and  Cj*  mode  coefficients 
demonstrate  a  sinusoidal  profile  with  the  notch  orientation  (Figures  6(a-b)). 


Figure  3:  Variation  of  fillet  radius  without  blind  damage:  (a)  &Sb  Mode  Coefficient 
(b)  Cjia  Mode  Coefficient 


Figure  4:  Variation  of  fillet  radius  with  blind  damage  of  constant  depth  and 
orientation:  (a)  $sa  Mode  Coefficient  (b)  Cta  Mode  Coefficient 


Figure  5:  Variation  of  blind  damage  depth  with  constant  fillet  radius:  (a)  ^sp  Mode 
Coefficient  (b)  Cj,q  Mode  Coefficient 


Figure  6:  Variation  of  blind  damage  orientation  with  constant  fillet  radius:  (a) 

Mode  Coefficient  (b)  ^Aa  Mode  Coefficient 

EXPERIMENTAL  RESULTS 

The  trends  from  the  FE  parameterization  are  investigated  by  using  a  3D  Scanning 
Laser  Vibrometry  test  setup  (Figure  7(a)),  which  consists  of  an  aluminum  T-joint 
for  selected  fillet  radius  (0.6  cm),  excitation  frequency  (100  kHz)  and  stiffener 
thickness  (0.3  cm).  The  setup  enables  3D  surface  measurement  of  the  front  and 
back  sides  of  the  plate,  plate-stiffener  radius,  and  stiffener  (Figure  (b)). 


(a)  (b) 


Figure  7:  Experimental  validation  using  Polytec  3D  Vibrometer:  (a)  Test  setup  (b) 
Aluminum  T-joint  specimen 

To  ensure  Lamb  wave  phenomena,  Figure  8(a)  illustrates  the  frequency  peaks  of  the 
incident  Ao  mode  of  the  front  face,  which  aligns  with  the  analytical  Ao  mode,  while 
Figure  8(b)  shows  the  back  face  frequency  content,  where  mode  conversion  occurs 
and  the  So  mode  is  introduced.  Figure  9(a)  demonstrates  the  full  wavefield  of  the 
filtered  out-of-plane  A0  mode  along  the  front  face,  and  Figure  9(b)  depicts  the  in- 


plane  converted  So  mode.  The  magnitude  of  the  fillet  radius  in  Figure  9(b)  is  3 
orders  below  that  of  the  front  and  stiffener  surfaces.  Using  data  whose  spatial 
points  originate  at  Y-axis  0.14  m  in  Figure  9(a)  and  runs  along  the  front,  back  and 
stiffener  faces,  Table  II  contains  the  reflection,  transmission,  and  conversion 
coefficients.  Particular  attention  is  given  to  the  difference  in  order  of  magnitude 
between  the  stiffener  and  front  and  back  faces,  indicating  that  a  significant  amount 
of  energy  is  propagating  into  the  stiffener.  The  relative  order  of  magnitude, 
approximately  0.02  to  0.06  along  the  front  face  for  the  reflection  and  conversion 
mode  coefficients  corroborate  with  the  results  from  the  undamaged  FE  trade  studies 
shown  in  Figure  3. 


Figure  8.  Incident  in-plane  frequency  content:  (a)  Front  face,  and  (b)  Back  face. 
Excitation  frequency  of  modulated  tone  burst  is  100  kHz. 


(b) 

Figure  9.  Time  history  snapshot  at  87.8  psec  along  front  face- stiffener  surfaces:  (a) 
Filtered  incident  out-of-plane  Aq  mode,  (b)  Filtered  converted  in-plane  So  mode 


Table  II.  Integrated  mode  coefficients  along  center  line  of  three  surfaces 


Section 

Reflection 

Coefficients 

Transmission 

Coefficients 

Conversion 

Coefficients 

D  U 

E-A0 

n  w 

Eao 

r>  uw 

Eao 

rp  U 

J-A0 

rp  W 

J-A0 

rp  UW 

J-A0 

pi  U 

t'SO 

pi  W 

t'SO 

pi  UW 

t'SO 

Front  Face 

0.054 

0.059 

0.058 

- 

- 

- 

0.053 

0.024 

0.039 

Back  Face 

- 

- 

- 

0.012 

0.019 

0.017 

0.014 

0.002 

0.009 

Stiffener 

- 

- 

- 

0.193 

0.203 

0.199 

0.118 

0.019 

0.074 

CONCLUSIONS 

As  the  fillet  radius  increases,  the  mode  coefficients  from  undamaged  to  damaged 
are  approximately  equivalent.  The  thickness  of  the  fillet  radius  (greater  than  6  mm) 
does  not  allow  the  damage  depth  to  penetrate  the  plate  thickness,  and  sensitivity  to 
the  relative  damage  is  therefore  negligible.  In  addition,  the  and  mode 
coefficients  demonstrate  a  sinusoidal  profile  with  the  notch  orientation,  such  that 
the  determining  the  orientation  based  on  the  magnitude  of  the  mode  coefficient  is 
limited.  Experimentally,  the  3D  Laser  vibrometry  results  indicate  that  the 
transmission  coefficients  are  an  order  of  magnitude  greater  within  the  stiffener 
section  than  the  front  and  back  faces.  Further  work  is  ongoing  in  varying  the 
damage  to  the  aluminum  T-joint,  and  results  will  be  presented  in  future  papers. 
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